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V  Abstract 

This  work  presents  an  examination  of  interpolative  schemes  for 
fatigue  crack  growth  at  elevated  temperature  (1200  F).  An  interpolative 
scheme  involving  the  linear  superposition  of  effects  due  to  sfross 
ratio,  loading  frequency,  and  peak  load  dwell  or  hold  time,  .is  applied 
to  two  state-of-the-art  craok  growth  rate  prediction  models.  These 
models  are  the  SINH  model,  developed  by  Pratt  and  Whitney  Aircraft,  and 
the  HSE  model,  developed  by  General  Electric.  SINE  is  based  on  the 

hyperbolic  sine  function,  and  MSE  is  based  on  a  sigmoidal  equation. 

C 

The  results  of  an  experimental  program  are  presented.  Fatigue 
crack  growth  rate  tests  are  performed  on  compact  tension  specimens 
according  to  ASTM  standards.  Two  interpolative  models  are  developed 
from  the  resulting  data  base.  Additional  tests  are  performed  and  the 
model  predictions  are  oompsred  to  these  additional  tests. 

It  is  found  that  linear  (on  log-log  scale)  functional  forms  can  be 

© 

assumed  to  relate  model  ooustants  to  the  test  parameters  (load  *atio, 
frequency,  and  hold  time).  Also,  it  is  found  that  these  functional 

^  forms  can  be  applied  to  both  models.  Uore  work  is  needed  to  determine 

whether  the  effects  of  variations  in  each  of  the  test  parameters  on 
craok  growth  rate  are  independent  of  each  other  and  can  be  linearly 
superimposed. 

When  the  same  functional  forms  are  applied  to  the  SINH  and  the  MSE 
models,  they  behave  much  alike. 

t. 
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EVALUATION  OF  INTERPOLATIVE  MODELING  CONCEPTS 
FOR  FATIGUE  CRACK  GROWTH  AT  ELEVATED  TEMPERATURE 


I.  Introduction 


Changes  have  been  reoently  made  in  the  United  States  Air  Foree 
design  and  maintenance  philosophies.  These  new  philosophies  have 
created  a  requirement  for  more  accurate  f at igue-oreck-growth-rate 
prediction  models  for  structural  alloys  used  in  modern  gas-turbine 
engines.  The  Engine  Structural  Integrity  Program  (ENSIP),  and  the 
Retirement  For  Cause  (RFC)  Program  are  two  Air  Force  Programs  which 
embody  these  new  philosophies  (1). 

ENSIP  specifier,  how  future  jot  engines  will  be  designed.  The  ENSIP 
design  philosophy  is  based  on  the  assumption  that  fatigue  cracks  exist 
in  new  engine  components  at  the  time  of  production.  It  then  becomes 
part  of  the  design  process  to  demonstrate  that  these  cracks  will  not 
grow  to  catastrophic  size  during  the  design  life  of  the  engine.  During 
the  1960's,  critical  engine  components  were  limited  by  creep  and  stress 
rupture  phenomena.  However,  the  lives  of  today's  critical  engine 
components  are  primarily  limited  by  low-cyole-fatigue  damage.  With  this 
new  emphasis  on  fatigue  phenomena,  crack  growth  prediction  has  become  an 
integral  part  of  the  engine  design  process. 

RFC  is  a  component-life  management  philosophy  which  seeks  to  safely 
utilize  a  greater  portion  of  the  lives  of  rotating  engine  parts,  and 
thereby  realize  considerable  economic  savings.  Consider  a  typical 
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turbine  engine  disk.  The  previous  approach  had  been  to  use  a  lower 
bound  on  the  Low-Cyc le-Fatigue  (LCF)  life  which  was  so  conservative  that 
it  resulted  in  999  out  of  1000  disks  being  retired  while  still 
structurally  sound.  Under  the  ourrent  RFC  approach,  disks  are  removed 
for  periodic  inspection  and  returned  to  service  if  no  cracks  larger  than 
a  specified  detectable  size  are  found.  Here  is  where  an  accurate 
oraok-growth-rate  prediction  capability  is  indispensable.  An  analysis 
must  be  able  to  demonstrate  that  cracks  of  the  detectable  size  and 
smaller  will  not  grow  to  catastrophic  size  during  the  next  inspection 
interval . 

Scone 

Most  current  fatigue  crack  growth  models  are  of  the  form 

da/dN  -  HAD  U) 
where  da/dN  is  the  crack  growth  rate  f.ud  4K  is  the  stress  intensity 
factor  range.  If  AK  is  to  be  a  valid  modeling  parameter.  Linear  Elastic 
Fracture  Mechanics  (LEFM)  must  hold  (2).  This  means  that  the  plastio 
zone  surrounding  the  orack  tip  must  be  small  when  compared  to  the  crack 
length.  At  elevated  temperatures  one  would  expect  time-dependent 
material  behavior  resulting  in  large  scale  plasticity  and  creep  damage. 
However,  several  investigators  (see  e.g.  Larsen  (1,3))  have  shown  that 
the  dK-aualysia  is  valid  for  nickel-base  superalloys  for  temperatures  up 
to  1350  F. 

The  goal  of  this  thesis  is  to  determine  the  effectiveness  of  state- 
of-the-art  models  in  predicting  crack  growth  rates  for  high  temperature 
applioations.  An  examination  of  two  representative  examples  of  current 
orack-growth-rate  predictive  models  which  use  4K  as  a  correlating  para- 
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meter  is  presented.  Both  of  the  models  are  ompirical  in  nature  since 
accuracy  requirements  for  design  exceed  the  predictive  capability  of  any 
existing  analytic  description.  The  two  models  examined  were  developed 
by  aircraft  engine  designers  Pratt  and  Whitney  (PW)  and  General  Electric 
(GE) . 


PW's  model  (SINH)  uses  the  hyperbolic  sine  function  to  fit  the 
basic  sigmoidal  shape  exhibited  by  da/dN  versus  4K  data,  while  GE's  uses 
a  modified  sigmoidal  equation  (MSE).  Both  models  are  discussed  in 
detail  In  Section  III, 

This  thesis  seeks  to  investigate  how  well  crack  growth  rate  models 
perform  when  they  have  been  developed  from  a  minimal  data  set.  The 
approaoh  taken  in  this  thesis  was  to  collect  a  minimal  set  of  data  at 
high  temperature  (120O4F),  varying  three  test  parameters:  stress  ratio, 
frequency,  and  hold  time.  Those  data  were  used  to  formulate  predictive 
models  for  both  the  SINH  and  the  MSE  equations.  Finally,  additional 
data  were  collected  and  compared  with  the  SINH  and  MSE  predictions  to 
assess  the  interpolativo  capabilities  of  the  models. 

It  has  been  demonstrated  that  both  the  SINH  and  the  MSE  equations 
do  a  reasonable  job  of  fitting  da/dN  versus  data  (4).  This  thesis 
differs  from  previous  work  in  that  it  attempts  to  test  the  validity  of 
models  in  which  more  than  one  of  the  test  parameters  are  allowed  to  vary 
at  any  one  time.  In  their  current  states  of  development,  both  the  SINH 
and  the  MSE  models  cannot  treat  these  oombined  effects  in  a  systematic 
manner  suitable  for  application  to  a  variety  of  materials.  The  version 
of  PW's  SINH  program  described  in  (5)  can  only  be  used  to  develop 
interpolat ive  models  in  whioh  only  one  test  parameter  is  varied  at  a 
time.  GE's  MSE  model  as  implemented  in  (6)  is  not  a  general  purpose 
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model.  While  it  allows  more  than  one  test  parameter  to  change  at  a 
time,  it  uses  ad  hoc  fractional  relationships  based  on  a  large  data  set 
for  one  specific  material.  Therefore,  it  does  not  attempt  to  produce 
general  functional  relationships  vhich  are  valid  for  a  range  of 
materials,  as  S1NH  does. 

In  this  thesis,  two  models  arc  developed  for  IN 718.  The 
experimental  program  used  in  the  development  of  these  two  models  is 
described  in  Section  II.  One  of  the  models  uses  the  SINH  equation,  and 
the  other  uses  the  MSB  equation.  Both  models  function  similarly.  Given 
a  set  of  test  parameters  (load  ratio,  R,  frequency,  V,  and  peak-load 
hold  time,  tjj)>  the  models  predict  da/dN  versus  4K  curves. 

Data  collection  for  fatigue  crack  growth  modeling  is  expensive  and 
can  take  several  days  to  perform.  Therefore,  in  order  to  reduce  the 
size  of  the  test  matrix  and  the  overall  complexity  of  the  modeling  task, 
the  temperature  was  held  constant.  All  data  were  collected  at  1200*F, 
Thus  temperature  was  eliminated  as  sn  independent  variable,  reducing  not 
only  the  total  number  of  test  runs  required,  but  also  the  number  of 
computer  runs  necessary  to  reduce  the  data  and  formulate  a  model. 

Fundamental  Assumptions 

In  both  the  SINH  and  the  USE  models,  da/dN  is  related  to  through 
a  set  of  constants;  SINH  uses  four  constants  while  MSE  uses  six.  When 
these  models  were  developed,  there  were  several  assumptions  made 
concerning  the  experimental  parameters  and  associated  constants.  These 
assumptions  are  suggested  In  the  SINH  docume ntat ion  (3),  but  are  not 
fully  implemented  in  the  SINH  computer  program.  They  have  never  been 
applied  to  the  MSE  equation  prior  to  this  work.  The  following 
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assumptions  are  made  in  this  thesis,  and  they  apply  to  both  the  S1NH  and 
the  MSE  models  for  IN718  developed: 

1.  Constants  are  related  to  the  test  parameters  through  rela¬ 
tionships  which  are  linear  when  plotted  on  log-log  graph 
paper. 

2.  Constants  are  functions  of  V,  tg,  E,  and  T,  where  V  is 
loading  frequency,  tH  is  load-dwell  time,  R  is  the  ratio 
of  minimum  to  maximum  load,  and  T  is  temperature. 

3.  First  order  dependencies  were  assumed,  i.e.,  it  is  assumed 
that  linear  superposition  holds,  and  that  there  are  no 
second  order  and  higher  effects. 

The  purpose  of  this  thesis  is  not  to  attempt  to  justify  these 
assumptions,  but  rather  to  examine  their  validity,  and  to  suggest 
possible  improvements. 


II.  Experimental  Program 


Computer-controlled  oonatant-load-amplitude,  fatigue  oraek-growth- 
rate  tests  were  performed  according  to  ASTM  Standard  E  647-81.  The 
experimental  program  was  performed  at  the  Air  Foroe  Materials 
Laboratory.  This  section  describes  the  experimental  apparatus, 
identifies  the  test  matrix  and  procedures,  and  outlines  the  data** 
reduction  techniques  employed. 

EMCilutmtil  Acputm 

A  closed  loop  computer-controlled  test  system  was  used  to  oolleot 
da/dN  versus  data.  A  block  diagram  of  the  test  setup  is  shown  in 
Figure  1.  The  primary  components  of  the  system  are  an  IBM  PC 
microcomputer,  an  MTS  Material  Test  Frame,  a  Wavatek  Function  Generator, 
and  various  Analog  to  Digital  (A/D)  converters  and  signal  smoothers. 
Data  are  taken  automatically  and  minimal  human  intervention  is  required. 

Figure  2  shows  the  two  raw  data  items  that  were  actually  measured 
by  the  test  set-up.  Crack  Opening  Displacement's  (COD)  are  measured 
using  an  MTS  electro-mechanical  extensometer  with  quart!  rods  for 
speoimen  contact,  while  the  loads  are  provided  by  an  MTS  servo-hydraulic 
load  frame,  as  controlled  by  the  computer.  The  specimens  are  enclosed 
in  an  electrical  resistance  furnace  to  achieve  the  required  elevated 
temperature . 

ISil  Matrix 

A  three-dimensional  representation  of  the  test  matrix  is  shown  in 
Figure  3.  Each  of  the  test  variables  is  plotted  on  a  separate  axis. 
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Figure  2.  Data  Acquisition  Prooess 
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Figure  3.  Graphic  Representation  of  Test  Matrix 
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The  four  points  plotted  as  solid  oiroles  represent  the  fonr  tests  used 
in  formulating  the  predictive  models.  The  coordinates  of  each  point 
are  the  values  of  the  test  parameters  for  that  particular  teat.  The  two 
points  plotted  as  squares  represent  the  two  proof  tests  which  were 
performed  to  evaluate  the  predictive  capability  of  the  models.  The  point 
at  the  origin  of  the  coordinate  axes  represents  the  baseline  test. 
Moving  from  the  baseline  outward  along  any  axis,  changes  only  the  value 
of  the  variable  represented  by  that  axis.  Thus,  the  three  points  lying 
on  the  axes  away  from  the  origin  represent  tests  in  whioh  only  one  test 
parameter  is  changed  from  the  baseline  value.  All  six  tests  were 
performed  at  a  temperature  of  1200  degrees  Fahrenheit. 

Specimen  Description 

Specimens  used  for  this  experimental  program  were  maohiued  from 
plates  of  INCONEL  718,  IN718  is  a  nickel-base  super  alloy  with  a  high 

fraoture  toughness,  and  is  commonly  used  in  the  fabrication  of  jet 
engine  hot-section  parts,  IN718  was  chosen  as  the  specimen  material 
beoause  of  its  utility  to  the  Air  Force,  and  beoauso  the  MSE  model  had 
previously  been  applied  only  to  AF115,  By  using  a  material  other  than 
AF115,  it  could  be  determined  if  the  functional  relationships  implicit 
in  the  MSE  model  could  be  applied  to  other  materials.  The  opecimen 
geometry  is  in  accordance  with  ASTM  Standard  E  647-81  (7),  and  is  shown 
in  Figure  4.  The  notoh  which  is  used  to  start  the  crack  was  cut  using 
Elootrioal-Discharge  Machining  (EDM). 

gmsdir.?. 

The  speoimens  were  pre-oracked  in  accordance  with  ASTM  Standard 
E  399  (8)  until  a  measurable  crack  extension  was  observed,  with 
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Figure  4.  Compect  Tension  (CT)  Specimen  Geometry 


the  final  pre-oraoking  Eaax  being  smaller  than  the  starting  Kmtx.  As 
suggested  in  ASTM  &  647-81,  Kmax  was  not  changed  by  sore  than  20%  in  any 
one  increment  so  as  to  establish  a  steady-state  crack  growth  rate  before 
taking  a  measurement.  Data  were  taken  automatically  with  the  equipment 
described.  Crack  length  was  derived  from  compliance  using  a 
relationship  obtained  from  a  polynomial  fit  to  a  seriea  of  numerical 
finite  element  solutions.  By  performing  a  series  of  finite  element 
calculations  for  the  compact  tension  (CT)  speoimen  at  different  oraok 
lengtha,  a  relationship  between  crack  length  and  compliance  oan  be 
obtained.  This  is  the  method  currently  in  use  at  the  Air  Force 
Materials  Laboratory.  For  a  discussion  of  the  error  inherent  in  this 
procedure  see  (9),  This  crack  length  expression  as  given  in  (9)  is: 

a  -  1.0010  -  4.66950"1  +  18.46U"2  -  236. 82U"3 

+  1214. 9D“4  -  2143. 6U“5  (2) 

where  U  -  (EBV/P)1/2  +  1  -  (EBC)1/2  +  1 
and 

a  -  average  crack  length  (inches) 

W  -  specimen  width  (inches) 

B  ■  speoimen  thickness  (inches) 

E  *  Young's  Modulus  (psi) 

V  -  crack  opening  displacement  (inches) 

P  “  load  ( inches /pounds) 

C  «  V/P  x  compliance  (pounds). 

The  craok  length-compliance  relationship  was  calibrated  at  the  beginning 
of  each  test  by  taking  a  few  visual  crack-length  measurements,  and  then 
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adjusting  Young's  modulus,  E,  in  Equation  2. 

The  stress  Intensity  range  4K  was  determined  from  the  ASTM 
expression  as  follows  (7): 


AV  2  +  a/W 
B(W)1/2  (1^  /W)3/ 


10.866  +  4.64(a/W) 


-  13.32U/W)2  -  5.62<a/W)3)  (3) 


where 

4P  -  load  range,  <PBax  -  Paia) 

Note  that  the  above  expression  assumes  a  linear-elastic,  homogeneous 
test  material. 

Pill  RtdmUga 

ASTM  Standard  E  647-81  suggests  two  methods  for  computing  crack- 
growth  rates:  the  Secant  method  and  the  Incremental  Polynomial  Method. 
The  Incremental  Polynomial  method  was  used  in  this  experimental  program, 
beoause  it  offers  some  smoothing  of  the  data.  In  this  method,  crack- 
growth  rates  are  computed  from  oraok  length  versus  number  of  cycles 
(a  versus  N)  data  using  a  seven-point  sliding  polynomial  ourve  fitting 
technique.  This  technique  gives  and  (da/dN)£^  by  the  following 
expressions: 

-  bQ  +  bjUNi  -  C^/Cj]  +  b2[(N1  -  C^/C^2  (4a) 

(da/dN)^  -  bj/Cj  +  2b2(Nl  -  C^/Cj2  (4b) 


’*  v-V 
**,  *  *  * 


, A", 

n.-'.-A 
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-1  iKNi  -  C^/CjU  +1 


«,  and  bj 


*» 


(da/dN)A  . 


oonstanti  daternined  by  least  squares  regression 
(Ni_3  +  Ni+3)/2 
(Ni+3  -  Ni_3)/2 
fitted  erack  length 

the  crack  growth  rate  aaaooiated  with  ^  and  N^. 
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m.  ffmth  IaUrngJ.fcU.Yc  MflAtliaa 

Typical  da/dN  versus  4K  data  plotted  on  log-log  graph  paper  hat  a 
sigmoidal  shape  aa  can  be  seen  In  Figure  5,  The  curve  ia  characterized 
by  a  alow  crack-growth  region  (I),  called  the  threahold  region,  a  linear 
region  (II).  and  a  rapid  growth  region  (III).  Data  do  not  nsually 
extend  into  region  III  becauae  teat  apeoinena  tend  to  break  aoon  after 
having  entered  thia  region. 

If  all  of  the  teat  paraaieters  were  held  constant  with  the  exception 
of  one.  aay  R,  variation  of  thia  parameter  would  generally  reault  in  a 
shifted  da/dN  veraua  dK  curve.  The  general  trends  for  these  curve 
shifts,  for  a  material  with  high  fracture  toughness,  such  as  IN718,  can 
ba  seen  in  Figure  6  (3,10).  Figure  6a  ahowa  that  increasing  v  will 
result  in  a  decrease  in  da/dN  for  a  fixed  value  of  AK,  Figure  6b  ahowa 
that  increasing  S  will  result  in  an  increase  in  da/dN  for  a  fixed  value 
of  /IK,  It  ia  the  goal  of  crack-growth  interpolative  models  to  be  able 
to  predict  these  ourve  shifts  with  reasonable  acouraoy,  from  a 
relatively  small  experimental  data  base.  What  follows  is  a  description 
of  the  two  interpolative  models  examined  in  thia  work,  and  how  they  were 
applied. 

SIM  MMcl 

Pratt  and  Whitney's  SINH  crack  growth  rate  interpolative  model  is 
based  on  the  hyperbolic  sine  function  (3,5),  as  plotted  in  Figure  7. 
Its  shape  closely  approximates  the  sigmoidal  shape  of  da/dN  versus  <lK 
data  shown  in  Figure  5.  With  the  addition  of  some  controlling  constants 
one  can  force  the  hyperbolic  sine  funotion  to  fit  may  possible 
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Crack  growth  rate  per  cycle,  da/dN  log  scale 
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variations  of  da/dN  versus  /\K  data.  The  SINH  modeling  equation  will  nov 
be  derived,  starting  with  the  definition  of  the  sinh  fnnetion. 


y  “  sinh  z  ■  - -  (5) 

2 

The  inflection  point  fox  sinh.  which  normally  falls  at  (0,0)  can  be 
shifted  by  the  addition  of  two  constanta  Cg  and  C4  as  follows: 

(y  ~  C4)  *■  s inh(-  +  Cg)  (6) 

Here,  Cg  locates  the  horizontal  position  and  C4  locates  the  ver'eioal 
position  of  the  inflection  point.  The  carve  can  be  stretched  vertically 
by  the  addition  of  a  constant  and  horizontally  by  the  addition  of  a 
constant  C2  such  that  the  equation  nov  reads 

(y  -  C4) 

-  -  sinh(C2(x  +  Cg),  (7) 

C1 

and  can  be  written  as 

y  *  C2ginh(C2<x  +  Cg)  +  C4  (8) 

Substituting  log(da/dN)  for  y,  and  log(^K)  for  z,  the  final  PW  SINH 
equation  becomes 

log  (da/dN)  -  CjsinhCCgdog^K)  +  Cg)  +  C4  (9) 

When  this  equation  is  applied  to  predict  crack  growth  rates,  the 
following  assumptions  are  made: 


V 


Q 


e 
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“  material  constant  (0.5  for  IN718) 

c2  “  f2<v'  tH*  R) 

^3  *  tH' 

C4  a  f4(  v>  tg,  H) . 

The  essence  of  the  SINH  model  is  to  find  these  functional 
relationships  which  relate  the  SINH  equation  constants  to  the  teat 
parameters.  In  this  work,  there  functional  relationships  were  all 
assumed  to  be  linear  on  log-log  plots.  Also,  it  was  assumed  that  the 
changes  due  to  v,  tg,  and  &  are  iudepeudent  of  each  other,  and  therefore 
their  effects  can  be  accounted  for  by  linear  superposition.  This 
suggests  that  the  values  of  the  constants  C2,  Cg,  and  C4  which  aooouut 
for  variations  in  each  of  the  test  parameters  are  obtained  from 

cj“  Cjba.eline+  4cj  J  -2.3.4  (10) 

where 

4C J  -  (dc./dv)Av  +  (dc.  /dtH)4tH  +  (dCj/6R)4R  (11) 

and  the  4's  refer  to  differences  from  the  baseline  values.  Specifically, 
whe2  j  •  2,  the  assumed  form  becomes 

h  ■  C2  baseline  H  V0**'  +  b2log(tH+l)  +  b3log( (1-R) / .9)  (12) 

where  the  last  three  terms  on  the  right-hand  side  acoount  for  the 
changes  in  C2  caused  hy  the  variation  of  each  test  parameter  from  the 
baseline  condition.  The  divisor  (.9)  in  the  last  term  forces  the  term 
to  vanish  at  the  baseline  value  of  R.  The  addition  of  unity  to  tg  in 
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the  next-to-last  term  is  a  slight  non-linearity  which  solves  the 
mathematical  difficulty  when  tg  equals  zero. 

The  constants  bj,  b2,  and  b3  appearing  in  the  last  three  terms  in 
Equation  12  were  determined  as  follows.  Referring  bach  to  Figure  3,  it 
can  be  seen  that  the  test  matrix  was  designed  to  allow  for  the 
segregation  of  the  V  .  tg,  and  K  dependencies.  There  are  two  different 
values  of  each  test  parameter.  This  forms  three  unique  groups  of  two 
data  sets.  Eaoh  of  these  groups  contains  the  baseline  data  set  and 
another  data  set  in  which  only  one  of  the  three  parameters  has  been 
varied.  Consider  the  group  of  data  sets  where  V  is  varied.  In  this 
case.  Equation  12  reduces  to  the  first  two  terms,  because  the  R  and  tg 
terms  vanish  at  the  baseline  values.  The  procedure  used  to  find  b^  was 
to:  (1)  fit  the  best  SINH  curve  to  the  baseline  data  set.  (2)  fit  the 
best  SINH  curve  to  the  set  in  which  V  is  varied,  and  (3)  fit  a  straight 
line  between  the  two  values  of  C2  found  in  steps  (1)  and  <2).  The 
slope  of  this  straight  line  is  bj.  The  same  process  was  repeated  for 
the  other  data-set  groups  to  obtain  the  last  two  terms  in  Equation  12. 
Expressions  for  the  remaining  SINH  constants.  Cj  and  C4.  were  obtained 
in  the  same  manner.  Pratt  and  Whitney’s  SINH  computer  program  was  used 
to  automate  this  process.  In  PW’s  SINH  computer  program  (5),  an 
iterative  soheme  using  least-squares  regression  is  used  to  find  the  SINH 
equation  constants  for  each  data  set  starting  with  initial  guesses 
supplied  by  the  user.  The  program  also  supplies  the  straight-lino  fits 
to  oonstants  derived  from  groups  of  data  sets  in  which  only  one  test 
parameter  is  varied. 
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MSB  Model 

The  Modified  Sigmoidal  Equation  (MSE)  model  was  developed  by 
General  Electric  (6),  and  is  given  by  the  following  expression: 

da/dN  -  eB'(4KMMi)P[lnC4i:Ml*)]Qtln0lKo/4K)]D  (13) 

Here 

■  At  at  the  inflection  point  of  the  ourve 
AK*  *  the  vertical  asymptote  at  crack  growth  threshold 

AK0  ■  the  vertical  asymptote  at  maximum  orsck  growth  rate 

Q  *  lower  shaping  coefficient 

D  »  upper  shaping  coefficient 

and 

D  -  -  [Q*^ln(4ljMK0) /ln(At^/At*)  ]2  (14) 

P  -  (da/dN1) '  -  d/lniAt^At*)  +  O/lniAtJAt^  (15) 

B'  «  ln(da/dNi)  -  ainUn^^/dS*)  3  -  Din  [In  (4^/4*^)  ]  (16) 

where 

da/dN ^  *  the  crack  growth  rate  at  the  inflection  point 

(da/dN^)'  ■  the  slope  at  the  inflection  point. 

Figure  8  shows  the  MSE  equation  plotted  and  indicate;  the  manner  in 
which  the  coefficients  in  the  MSE  equation  affect  the  la/dN  versus  At 
curve  ( 6 ).  The  coefficients  Q  and  D  control  the  sharpness  of  the 
transition  from  the  inflection  point  to  the  asymptotes.  Decreasing  the 
absolute  value  of  these  constants,  causes  a  sharper  transition.  When  D 
is  set  equal  to  minus  Q,  a  symmetric  curve  results.  B*  is  the  distanoe 
by  which  the  inflection  point  is  displaced  from  a  da/dN  value  of  one. 
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Figure  8.  The  Effect  of  the  MSE  Coefficients  os  the  Siguoidsl  Curve 

(Frou  GE  Report  (6)) 
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The  unprimed  B  in  the  figure  oontrols  the  vertical  motion  of  the  whole 
carve.  P  controle  the  rotation  of  the  carve  in  the  vicinity  of  the 
inflection  point. 

A  comptriion  of  the  SINH  end  MSB  models  reveals  that  there  are 
analogies  between  the  two  sets  of  coefficients  employed  by  each  model. 
Coefficients  da/dN^,  (da/dN^)',  4K0*  At* ,  and  Q  are  analogous  to 
the  constants  Cj-C4  in  the  SINH  model.  Koto  that  the  constants*  At±  and 
da/dN^,  are  related  directly  to  the  two  SINH  constants*  Cg  and  C4  by  the 
following: 


C3  -  -log^Kj)  (17) 

and 

C4  -  log(da/dNl)  (18) 

As  it  was  done  with  the  SINH  model*  expressions  were  derived 
relating  these  six  constants  to  the  test  parameters.  The  procedure  used 
to  determine  the  SINH  constants  was  also  followed  in  determining  the  MSE 
constants.  However*  since  GE's  MSE  model  is  not  developed  to  the  point 
of  a  production  code,  any  regression  analysis  needed  was  performed  by 
hand  calculations.  In  fact,  as  mentioned  earlier,  the  GE  MSE  model  was 
specifically  developed  for  the  alloy  AF115,  and  has  not  been  previously 
applied  to  IN718,  The  General  Eleotrio  investigators  (6)  had  aooess  to 
a  large  data  set,  and  thereby  were  able  to  include  some  higher  order 
dependencies,  i.e.  terms  involving  products  of  functions  of  the  test 
parameters.  The  data  set  used  in  this  investigation  was  limited  by 
dosign.  Therefore,  the  linear  forms  used  for  the  SINH  model  were  also 
applied  to  the  MSE  model.  One  of  the  goals  of  this  investigation  is  to 
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determine  whether  the  nee  of  linear  functional  forme  ie  adequate.  The 
assumed  forme  are  given  as  follows: 

4K0  -  d2  (19) 

logOlK*)  -  log(4K*)^t8e  +  djlogl/  +  d j log ( tg+1 )  (20) 

logOlX.^)  ■»  logOlX^)^,,  +  d^log(l-R)  +  djlogt'  +  dglog(tH+l)  (21) 

log(da/dN^)  ■  log(da/dN£)^tte  +  d<ylog(l~R)  +  dglogP  +  dplog(tjj+l)  (22) 

<da/dNi)'  -  (da/dNj)  'ba$e  +  d10log(l-R)  +  dj^logp  +  d^ogU^l)  (23) 

Q  -  d13  (24) 

where  dj  through  dj3  are  constants  determined  by  linear  regression, 
while  is  determined  graphically  via  an  iterative  prooees.  There  are 
only  two  data  curves  need  to  resolve  a  change  due  to  varying  a 
particular  test  parameter  in  Equations  19-24.  Therefore,  the  linear 
regression  approach  reduces  to  solving  a  set  of  simultaneous  equations. 
For  example,  in  Equation  22,  the  set  of  simultaneous  equations  is  a 
three  by  three  set.  The  value  of  da/dN^  is  known  for  the  baseline  data 
set  and  for  the  three  data  sets  where  V,  R,  and  tjj  are  varied.  These 
values  can  be  substituted  into  Equation  22  along  with  the  appropriate 
values  of  the  test  parameters  for  each  data  set.  The  result  is  three 
equations  in  the  three  unknowns  d^ ,  dg,  and  dp. 

The  following  comments  should  be  made  concerning  these  assumed 
functional  forms:  (1)  The  upper  asymptote,  AKc ,  was  assumed  to  be  a 
constant  for  simplicity  and  because  the  experimental  data  did  not  extend 
into  this  region.  This  does  not  over-restrict  the  shape  of  the  curve, 
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because  if  Q  i»  not  equal  to  minn#  D,  the  bottom  half  can  act 
independently  of  the  top  half.  (2)  An  examination  of  the 
experimental  data  revealed  that  AK*  ia  only  a  function  of  v  and  tfl  and 
not  R.  (3)  The  elope  of  the  curve  at  the  inflection  point,  (da/dN^', 
vae  obtained  visually  from  data  plots.  The  points  determined  using  the 
SINH  computer  program  to  be  the  inflection  points  for  each  data  set, 
vere  also  used  for  the  MSE  model.  Therefore,  the  inflection  points  of 
curves  produced  from  both  models  coincide.  (4)  The  lower  shaping 
coefficient,  Q,  was  assumed  to  be  a  constant,  since  it  was  found 
graphically  via  an  iterative  process  that  adequate  MSE  fits  oould  be 
obtained  for  all  of  the  experimental  data  sets  with  a  oonstant  Q.  The 
first  attempt  to  model  Q  required  the  use  of  an  expression  similar  to 
that  given  in  Equation  23.  However,  this  resulted  in  an  expression 
which  gave  unreasonable  values  of  Q  for  values  of  the  test  parameters 
corresponding  to  the  proof  teats,  i«e„  negative  Q's.  It  was  oonclcded 
that  Q  was  sensitive  to  change  near  the  baseline,  but  was  relatively 
insensitive  to  values  of  the  test  parameters  away  from  the  baseline. 
Therefore,  the  value  of  Q  fo’  ed  to  best  represent  the  baseline  condition 
was  also  used  for  the  remaining  data  sets.  (5)  The  motivation  for 
assuming  MSE  equations  of  a  form  analogous  to  that  used  for  the  SINH 
model,  was  to  determine  whether  they  are  valid  for  either  model.  They 
are  general  equations  which  oould  be  applied  to  either  a  small  or  a 
large  data  set.  The  objective  was  to  produoe  functional  relationships 
which  do  not  rely  on  ad  boo  relations  suggested  by  the  data  itsolf,  but 
remain  valid  for  any  material  exhibiting  similar  fracture  toughness 
characteristics. 
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IV.  Results  ifli  Discussion 


Experimental  Bill 

All  of  the  data  obtained  from  the  experimental  program  are  given  in 
numerical  form  in  Appendix  A.  The  data  are  grouped  aooordiag  to 
specimen  numbers  assigned  at  the  time  of  testing.  The  four  data  sets 
used  for  developing  the  SINH  and  MSE  models  are  given  first.  Their 
specimen  numbers  are  832(540  (baseline  data  set,  R-.1.  v -1  Ha,  tg-0  sec), 
83266G  (R-.l,  •'■•01  Hx,  tH-0  seo),  83270B  (R-.5,  ^-l  Hr,  tH-0  seo),  and 
83270  (R-.l,  v-1  Hx,  tH-S0  seo).  The  following  two  data  sets  were 
obtained  from  the  proof  tests.  These  specimens  were  used  to  check  the 
predictions  of  the  models  developed  from  the  first  four  data  sets.  The 
specimen  numbers  are  83288  (first  proof  test,  R—.5,  *'—•01  Hx,  tjj-50 
seo),  and  83289  (second  proof  test.  R-.25,  v-1  Hx,  tg-10  sec).  Finally, 
there  are  two  groups  of  constant  48  data  provided  from  the  archives  of 
the  Air  Force  Materials  Laboratory  taken  from  tests  performed  prior  to 
this  thesis  effort.  The  first  of  these  sets  is  for  R-.l,  tH**0  sec,  and 
various  frequencies  at  48-22.745  Esi  Vinches  and  48-36.4  KsiVinoheT. 
The  seoond  of  these  two  sets  is  for  R-.l,  v-1  Hx,  and  various  hold  times 
at  the  same  two  4K's.  These  two  data  sets  were  used  to  examine  the 
validity  of  the  linear  functional  relationships  assumed  for  v  and  tg. 

The  data  used  for  the  development  of  the  SINH  and  MSE  models  are 
plotted  in  Figures  9-12  and  also  in  Figures  13-16.  In  Figures  9-12,  the 
dashed  lines  are  the  SINH-modol  fits  to  the  data  sets,  while  in  Figures 
13-16  the  dashed  lines  represent  the  MSE  fits  to  the  data.  The  large 
solid  squares  in  Figures  9-12  identify  the  locations  of  the  inflection 
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point*  calculated  iteratively  by  PW's  SINH  program  atarting  with,  uaor 
supplied  guesses.  It  should  be  noted  that  tho  final  looatlon  of  the 
inflection  point  depends  upon  the  user's  initial  guess.  This  point  is 
discussed  in  detail  in  (4)  . 

SINH  and  MSB  Predictions 

The  reader  will  recall  the  prooodure  described  in  Section  III 
regarding  the  segregation  and  characterization  of  the  changes  in  t.\e 
SINH  (or  NSE)  constants  due  to  changes  in  the  test  parameters  R,  v,  and 
tg.  This  procedure  consists  of  grouping  experimental  data  sets  in  which 
the  values  of  all  the  test  parameters  are  the  sane  except  for  one. 
Equations  are  then  developed  which  allow  for  linear  interpolation 
between  the  model  constants  obtained  for  each  data  set  in  a  group. 
After  the  effects  ou  the  constant*  due  to  changes  in  esoh  of  the  test 
parameters  ara  isolated  and  characterised,  they  are  superimposed  to  form 
e  single  equation  for  eaah  model  constant.  Figure*  17-19  illustrate  how 
this  procedure  was  executed  for  the  SINH  model.  Figure  17  shows  the 
change  in  the  position  of  the  inflection  point  due  to  a  change  in  R, 
while  Figures  18  and  19  show  the  changes  due  to  changes  in  v  and  tg 
respectively.  In  these  figures,  the  solid  lines  drawn  through  the 
inflection  points  represent  the  sssooiated  linear  fits.  The  dashed 
lines  ere  the  model  fits  to  each  date  set.  In  each  of  the  Figures  17- 
19,  only  one  teat  parameter  is  varied.  Similar  lines  could  be  drawn  for 
the  MSE  fits.  In  fact,  they  would  be  the  same  lines  beoause  the  same 
inflection  points  wore  ohosen  for  the  SINH  and  MSE  curves  for  each  data 
aet,  even  though  the  shape  of  the  curves  is  slightly  different. 

Using  tho  procedures  outlined  in  Section  III,  the  following 


e 


expressions  were  derived  reproseating  appropriate  interpolativo  models 
for  the  SINH  sad  the  MSB  equations: 

SINH  Intoroolative  Mflifil 

log.da/dN)  *  C^einh(C2( log(4K)  +  C3)  +  (25) 

where 


M2E  Interoolatlve  Model 

da/dN  -  eB'(4KM«:i)P[ln(4K/aK:*)]Q[ln(4K0MK)]D 


whe^e 


4K0  -  109.9 


4K*  -  iq[1. 08918  -  ,048461ogt'  +  .056751og(te+l) ]  (32) 

AK£  -  io II *721  +  . 6243 log ( (1-R) / . 9)  -  ,007171oge  -  .O6771og( tH+l> ]  (33) 
Q  M  1.0  (34) 


da/dNj.  -  10 


[-3.935  +  1. 0131 log ( (1-R) /.9) 


■  .5 

(26) 

M 

-  3.197  -  1.1675 log [(1-R)/. 9]  - 

.15051ogi/  +  .20531og(tH+D 

(27) 

‘(.  «\\Vj 

-  -1.721  -  .6243 log [(1-R) / .9]  + 

.00717  log v  +  .06771og( tg+1) 

(28) 

O 

■jv  jv.% 

-  -3.935  +  1.0131 log [(1-R0/ .9] 

-  .523  log v  +  . 8127 log (tH+l) 

(29) 

-  .  5231og  v  +  .  8127  log  ( trj+1)  ] 
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D  "  -  [Q1/2ln(4Ki/AK0)/ln<4KiMK*)]2  (36) 

P  -  (da/d^)'  -  Q/ln(4Ki/AK*)  +  D/  la(AK0MK^)  (37) 

B'  -  InUa/dN^  -  QlntlnCdK^dK*) )  -  Dln[ln(4K(J/4Ki)] .  (38) 

(da/dNi)*  “  1.7  -  1.84511og[(l-&)/.9] 

-  .21851og  +  , 1956 log ( tg+1)  (39) 


Note  that  Equation*  25  and  30  were  previously  given  as  Equationa  9  and 
13  respectively.  Also,  Equations  36,  37,  and  38  are  the  same  as  14,  15, 
and  16.  In  the  above  expressions,  C3  has  the  units  of  log(Ksi  VfiT.),  C4 
has  the  units  of  log( inches/ cyo le) ,  C^,  C2,  (da/dN^)',  Q,  D,  and  P  are 
dimensionless,  4K0.  Ah i ,  and  dK*  have  the  units  of  Ksi  vTn*.,  and  B'  has 
the  units  of  lnUnohes/cycle).  Both  the  SINE  and  the  MSE  Interpolative 
Models  were  programed  in  standard  FORTRAN,  and  the  resulting  interactive 
computer  codes  appear  in  Appendix  B  and  in  Appendix  C,  respectively. 
These  similar  codes  calculate  the  constants  for  the  corresponding  model 
at  T»  1200  F,  when  R,  V,  and  tg  are  given  as  input.  They  also  give 
da/dN  values  for  requested  dK  values.  They  are  based  on  a  limited  data 
set,  and  they  are  given  here  to  show  how  the  models  wore  applied. 

Pl89ttAAl9fl 

A  few  comments  are  in  order  concerning  the  SINH  and  MSE 
interpolative  models  as  developed. 

In  Equations  27,  28,  29,  33,  35,  and  39,  the  .9  divisor  in  the  R- 
ohange  terms  has  been  included  as  an  offset  used  to  force  the  R-ohange 
term  to  vanish  when  R  is  at  its  baseline  value. 
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Figure*  13-16  show  the  initial  USE  fits  to  the  four  data  sets  from 
which  the  model*  were  constructed.  These  fits  were  obtained  by  hand 
calculation*.  However,  the  initial  SINH  fits,  shown  in  Figures  9-12, 
were  obtained  with  the  sid  of  Pratt  and  Whitney's  SINH  computer  program 
which  uses  an  iterative  regression  scheme.  Even  with  this  aid,  a  great 
deal  of  user  judgement  is  required  in  order  to  obtain  a  good  fit.  The 
initial  SINH  fits  were  obtained  first.  Then,  the  USE  constants  which 
control  the  inflection  point  were  given  analogous  values  found  from  the 
initial  SINH  fits.  With  the  inflection  point  identified  in  the  MSE 
model,  and  using  a  oonstant  lower-shaping  coefficient,  Q,  and  assuming 
that  the  two  asymptotes  are  correct,  the  only  quantity  that  can  improve 
the  initial  USE  fits  is  the  oonstant  (da/dNp*.  Beoall  that  this 
oonstant  controls  the  slope  of  the  curve  at  the  inflection  point.  All 
of  these  initial  USE  slope  constants  were  computed  by  hand,  by  visually 
measuring  the  slopes  from  data  plots.  In  order  to  get  the  best  possible 
USE  fits,  it  was  necessary  to  adjust  both  Q  and  (da/dNp’ 
simultaneously.  Using  a  larger  slope  at  the  inflection  point  requires 
the  use  of  a  smaller  Q.  in  order  to  keep  the  same  ourvature  in  the  lower 
portion  of  the  curve. 

The  interpolative  models  were  compared  against  two  proof  tests. 
The  data  from  these  proof  tests  were  not  used  in  the  formulation  of  the 
models,  and  were  taken  at  arbitrary  values  of  the  test  parameters.  The 
results  obtained  from  Proof  Test  One  (R“.5,  I'-.Ol  Ha,  tg**50  sec)  are 
shown  in  Figure  20,  while  those  from  Proof  Test  Two  (R">,25,  v  -1  Hz, 
tgw10  toe)  are  shown  in  Figure  21.  Let  us  first  examine  Proof  Test  One 
(Figure  20).  Both  the  SINH  and  the  USE  models  matched  closely  in  their 
predictions,  with  the  MSB  predicting  s  steeper  slope.  SIND  appears  to 
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Figure  21,  Comparison  of  Predicted  and  Experimental  da/dN 
Proof  Teat.  Two  (tt-,25,  v-l  H* ,  tg-10  sec) 


t' 


m 


o 


0 


c 


predict  a  more  oorreot  threshold  4k  (SINH  does  not  have  a  true 
threshold)*  hut  both  seemed  to  predict  da/dN  values  greater  by  a  factor 
of  about  four  than  the  experimental  results.  Also,  the  slope  of  both 
predicted  curves  seemed  to  be  too  steep. 

Consider  next  the  results  of  the  Second  Proof  Test  (Figure  21). 
Here,  the  two  predictions  match  each  other  closely  except  for  the 
natural  differences  between  the  two  equations.  Also,  the  threshold 
value  of  4K  was  predicted  more  closely  than  in  the  First  Proof  Test. 
The  overall  aoouraoy  of  the  prediction  was  much  closer,  say.  within  a 
factor  of  one  to  two.  This  better  aoouraoy  may  be  attributed  to  the 
faot  that  the  values  of  the  teat  parameters  for  these  data  are  olosor  to 
the  baseline  values.  Also,  only  two  test  parameters  were  changed  from 
the  baaeline  values,  while  all  three  wore  changed  in  Proof  Test  One, 
The  limited  accuracy  of  these  two  predictive  comparisons  is  not 
surprising  considering  the  variability  Inherent  in  fatigue  oraek-growth 
data.  ASTM  standard  E  647-81  (7)  states  that  da/dN  data  may  vary  by  a 
factor  of  two.  Considering  that  only  two  data  points  were  used  to 
establish  each  linear  relationship  in  a  system  given  to  those  kind  of 
variubilities,  it  is  not  surprising  that  the  predictions  could  be  in 
error  by  a  factor  of  four. 

Another  type  of  oompnrison  was  made  using  the  constant  4K  data 
given  in  the  last  two  pages  of  Appendix  A.  Figures  22  and  23  compare 
da/dN  for  various  hold  times  at  4K>32.75  Ksi  VTn.  and  41022.724  Ksi 
VTn".,  respeotively.  This  data  was  taken  from  a  different  batch  of 
specimens,  whioh  explains  the  offset  of  approximately  a  faotor  of  two. 
Figures  24  and  25  compare  da/dN  for  various  frequencies  also  at  41032.75 
KsiVTn.  and  41022.745  Ksi  VFiT.  .  Figures  22-25  represent  an  attempt  to 
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Figure  22.  CoupMi'ieon  of  Predicted  tad  Experimental  da/dN  Date 
for^I-32.75  Koi  root  in  (H-.l,  ^-1  «*»  Variable  tH) 
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Figure  24.  Comparison  of  Predicted  and  Experimental  da/dN  Data 
for  ill-32.75  Ksi  root  in  (R“.lf  tg»0  see,  variable v  ) 
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FREQUENCY  (HERTZ) 


Figure  25,  Comparison  of  Predicted  and  Experimental  da/dN  Data 

for  ^£•‘22.745  Kai  root  in  (R-.l.  tg“G  sec,  variable  V) 


examine  the  validity  of  assuming  linear  functional  relationship*  (on 
log-log  plots)  relating  the  test  parameters  to  the  SINH  and  MSE 
constants.  In  these  plots  only  one  test  parameter  vas  varied  at  a  time, 
with  the  other  test  parameters  kept  at  baseline  values.  Figures  22  and 
23  show  da/dN  versus  hold  time  (tg)  for  two  different  4K's.  Notice  that 
the  SINH  and  MSE  curves  aro  slightly  non-1 inoar.  Aa  mentioned  earlier, 
this  is  because  log(tg+l)  was  assumed  as  a  functional  relationship 
rather  than  the  log(tg)  as  was  suggested  in  C3).  Assuming  this  form, 
takes  care  of  the  mathematical  difficulty  encountered  at  tg«*0,  and  seems 
to  match  the  form  of  the  experimental  data  better  than  a  straight  line 
would.  This  is  true  at  least  for  IN718.  In  Figures  25  and  26,  da/dN 
versus  frequency  (V)  is  plotted  for  two  different  <dK*s.  Note  that  in 
both  figures  the  experimental  data  shows  a  definite  upward  turn  at  the 
lower  frequencies.  This  cannot  be  modeled  by  a  straight  line.  However, 
considering  the  amount  of  variability  of  the  experimental  data,  a 
straight  line  seems  adequate  for  a  wide  range  of  V. 


V.  Mid.  Recommendations 


Form  al  Craok  Growth  Site  iaMlifitt 

Thin  exercise  ha*  shown  that  it  make*  little  difference  which  fora 
of  equation  io  used  to  model  fatigue  crack  growth  rate  data.  Both  the 
Hyperbolic  Sine  and  the  Modified  Sigmoidal  Equation  have  sufficient 
flexibility  to  model  da/dN  data.  Both  forms  give  essentially  the  same 
answers  when  the  same  functional  forms  are  used  to  relate  the  test 
parameters  to  the  equation  constants.  The  MSE  equation  offers  a  little 
more  flexibility  in  that  a  non-symmetrio  curve  can  be  modeled  as  well  a* 
a  symmetric  curve.  Also,  the  MSE  equation  has  the  advantage  of 
exhibiting  true  asymptotic  behavior  at  the  threshold  and  high  oraok- 
growth-rate  regions  of  the  da/dN  versus  curve. 

Xt  has  been  shown  that  the  same  functional  relationships  used  in 
the  SXNH  model  oan  be  applied  to  the  MSE  model.  Therefore,  the  same 
iterative  regression  techniques  employed  in  Pratt  and  Whitney's  SINH 
computer  program  could  be  applied  to  the  MSE  model.  Thus,  the  MSE  model 
could  be  automated,  making  it  a  more  usable  tool. 

Relating  Constants  i a  Test  ggxamtAXJ. 

The  method  used  to  relate  the  SINH  or  MSE  constants  to  the  test 
parameters  is  what  makes  either  one  of  them  a  model  and  not  just  an 
equation.  Certainly,  if  every  possible  combination  of  test  parameters 
were  tested,  then  relationships  could  be  developed  to  model  exactly  auy 
coudition  of  interest.  However,  if  all  possible  tests  were  performed, 
tbe  model  would  no  longer  be  needed.  The  reason  for  the  model  it  to 
reduoe  the  number  of  tests  necessary  to  characterise  the  oraok  growth 


behavior  of  the  aaterial  in  question.  Even  if  all  possible  testa  were 
performed,  there  would  still  be  the  possibility  of  a  two  hundred  peroent 
error  from  test  to  test  performed  under  the  same  conditions.  The 
question  of  whether  linear  (on  log-log  plots)  relationships  are  adequate 
to  relate  the  model  constants  to  the  test  parameters  cannot  be  fully 
explored  with  such  a  small  data  aet  as  was  used  here.  However,  any 
higher  order  interactions  would  probably  yield  changes  in  da/dN  which 
are  much  smaller  than  the  200  peroent  variability  of  the  data. 
Therefore,  the  question  is  probably  moot.  If  linear  (on  log-log  plots) 
relationships  are  employed  to  oonatruot  SINH  or  HSE  models,  the  models 
are  adequate  for  a  large  range  of  diX,  This  is  true,  providing  the 
region  of  model  use  is  not  too  far  from  the  baseline  condition.  The 
reader  will  reoall  that  it  was  assumed  that  changes  in  da/dN  due  to 
variations  in  eaoh  of  the  test  parameters  are  independent  and  can  be 
linearly  superimposed,  The  model  would  be  mo, re  accurate  near  the 
baseline  condition  because  any  inacouraoies  introduced  by  this  linear 
superposition  would  be  minimized  when  summing  small  individual  changes. 
A  larger  number  of  tests,  with  replication  of  as  many  tests  as  possible 
is  desirable  in  order  to  develop  any  real  confidence  in  models  developed 
from  a  particular  data  base. 

Future  work  should  investigate  the  type  of  functional  relationships 
that  are  appropriate,  using  a  much  larger  data  base  than  tbo  one 
generated  for  this  investigation.  Based  on  the  examination  of  the 
results  obtainod  from  the  two  proof  tests,  it  would  seem  that  linear 
superposition  of  slopes  and  inf leo tion-point  coordinates  might  not  be 
valid.  On  the  other  hand,  with  such  a  small  data  base,  small  errors  in 
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Specimen 

83264D  (Bueline  Pits 

Set.  R*.l»  V 

di/dN 

AIL 

12.599000 

0.61370000E-06 

22.735200 

12.625400 

0. 14596000E-05 

23.097300 

12.678200 

0.20198000E-05 

23.469500 

12.758300 

0.24883000E-05 

23.890700 

12.832000 

0.24217000E-05 

24.347400 

12.885600 

0.20665000E-05 

24.855100 

12.927500 

0.23672000E-05 

25.391000 

12.992100 

0.25621000E-05 

25.867700 

13.064000 

0.27196000E-05 

26,373600 

13.154000 

0.26661000E-05 

26.938600 

13.226800 

0.29269000E-05 

27.534500 

13.305100 

0.30447000EH95 

28.200500 

13.372400 

0.30093000E-05 

28.924700 

13.467000 

0.33797000E-05 

29.758100 

13.582600 

0.32559000E-05 

30.739800 

13.680800 

0.34933000E-05 

32.005300 

13.764500 

0.32903000E-05 

32.076300 

13.852800 

0.36374000B-05 

32.815000 

13.966500 

0.37328000E-05 

33.522900 

14.093000 

0.37527000E-05 

34.207900 

14.225800 

0.38424000B-05 

34.964900 

14.325900 

0.34922000E-05 

35.816500 

14.418700 

0.38277000B-05 

36.762700 

14.326000 

0.39996000E-05 

37.800800 

14.676100 

0.44189000E-05 

38.973500 

15.000900 

0.48992000E-03 

40.311900 

15.159300 

0, 52055000 E-05 

41.823100 

15.324800 

0 . 52990000E-05 

43.530800 

15.516800 

0. 57335000E-05 

45.486000 

15.734300 

0.61535000E-05 

47.642200 

15.965300 

0.63996000E-05 

16.207400 

0.68929000E-05 

16.470300 

0.70933000E-05 

16.766000 

0.75638000E-05 

17.071700 

0.79354000E-05 

17.407400 

0.85925000E-05 

18,186200 

0.96055000E-05 

18.639300 

0.10954700E-04 

19.184300 

0.12741300E-04 

19.858400 

0.15042900E-04 

20.732800 

0.18056700E-04 

20.774600 

0.15779500E-04 

21.016600 

0.15861000E-04 

21.220400 

0.15572400E-04 

21.457900 

0.16617300E-04 

21.745400 

0.18522000E-04 

22.053800 

0.19684600E-04 

,  tjj-0  teo) 
da/dN 

0.21055500E-04 

0.21737800E-04 

0.22763300E-04 

0.23882600E-04 

0.25173600E-04 

0;26543600E-04 

0.28223600E-04 

0.294833008-04 

0.30638900E-04 

0.32263300E-Q4 

0.34285000E-04 

0.36271200H-04 

0.38646000E-04 

0.44244000E-04 

0.51129800E-04 

0.56401500B-04 

0.54173100E-04 

0.55940800E-04 

0.55535300E-04 

0.536B4900E-04 

0.57184900E-04 

0.61267600E-04 

0.65133700E-04 

0. 89623 800E-04 

0. 74743 900E-04 

0.79684900&-04 

0.8467 7000 E-04 

0.90114000E-04 

0 95054900E-04 

0.10019670E-03 
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Specimen  832 660 
da/dN 

13.793400  0.55646000E-05 
15.804300  0.16774000E-04 
13.826100  0.27383000E-04 
15.838000  0.38392400E-04 
13.918000  0.47370000E-04 
13.976300  0.46582600E-04 
16.000800  0 , 28948800E-04 
16.026300  0.33642100E-04 
16,060000  0.40909000E-04 
16.143700  0.50078600E-Q4 
16.212800  C, 48873900 E-04 
16.272900  0.44527500E-04 
16.298300  0.39444800E-04 
16.327500  0.39862100E-04 
16.365700  0.38914500E-04 
16.413900  0.41507800E-04 
16.465700  0.48637700E-04 
16.513100  0.51602300E-04 
16.567600  0.40893600E-04 
16.609500  0.40783400E-04 
16.649500  0.47822700E-04 
16.701400  0.45873900E-04 
16.747800  0.39123200B-04 
16.785100  0.43031400E-04 
16.848800  0.46992000E-04 
16.877000  0.30397500E-04 
17.013400  0.54121900 E-04 
17.079900  0.53224300E-04 
17.121700  0.5376534OE-O3 
17.226300  0.559054002-04 
17.302800  0.65523500E-04 
17.403800  0.70507700E-04 
17.479300  0.62114000E-04 
17.536600  0.56098300E-04 
17.600300  0.52549000E-04 
17.700400  0. 53602200E-04 
17.775000  0.53066800E-04 
17.830300  0.60535300E-04 
17.393200  0.72212500E-04 
17.997900  0.81842400E-04 
18.107900  0.82322700E-04 
18.215300  0.79275400E-04 
18.326300  0.80759700E-04 
18.390000  0 , 80094300E-04 
18.479200  0.81177000E-04 
18.596500  0. 834565 00E-04 
18.715700  0.80614000E-04 
18.796700  0.82696700E-04 
18.867600  0 . 83161300E-04 
18.921300  0.84527400E-04 


V».01  Hi,  tg-0  sec) 

4K  da/dN 

19.133300  0.10171240E-03 

19.239800  0.99380000E-04 

19.414400  0.87531300E-04 

19.488100  0.89881700E-04 

19.580900  0.10008640E-03 

19.687400  0.10799580E-03 

19.817500  0 . 11166120E-03 

19.948500  0.11499980E-03 

20.056800  0.10773600E-03 

20.165000  0.10433440E-03 

20.268700  0.10581470B-03 

20.375200  0.11861790E-03 

20.518900  0.12694860E-03 

20.669100  0.13047000K-03 

20.830100  0.13185010E-03 

20.981100  0.13275560E-03 

21.121200  0.13344070E-03 

21.265900  0.13142490E-03 

21.436000  0.13746430E-03 

21.602500  0 . 14221800E-03 

21.777200  0.14862180E-03 

21.950100  0.15641310E-03 

22.153000  0.15969000E-03 

22.364900  0.17043270E-03 

22.602400  0.16697210E-03 

22.819900  0.16228000E-03 

23.161000  0.15706270E-03 

23.305700  0.16149970E-03 

23.424900  0.17466890E-03 

23.608700  0.20541300E-03 

23.818800  0.21221220E-03 

24.049900  0.22533030E-03 

24.266500  0.21643260E-03 

24.489400  0.21513340K-03 

24.672200  0.21090510E-03 

24.884200  0 .21170040E-03 

25.109800  0.22623940E-03 

25.345500  0.23283420E-03 

25.603900  0.24760190B-03 

25.874100  0.26016480E-03 

26.164300  0.26456250E-03 

26.481800  0. 27666480E-03 

26.795700  0.27952700E-03 

27.127800  0.28825930E-03 

27.470800  0.29139310E-03 

27.818400  0.30559780E-03 

28.629900  0.33940880E-03 

29.099400  0 . 34902690E-03 

29.584300  0.36035360E-03 

30.091100  0.37663310E-03 
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• 

Specimen  83270B  (R-.5,»/-l  Hi,  tg-0  sec) 

dt/dN  £&  <*a/dN 

“  .  '  '  •  ‘.4 
*..**  .‘*.1 

.  ■  *,y 

26.770200 

0.33575100E-04  12.079500  0.51524000E-05 

23.921600 

0 .28068100E-04  12.210600  0.54268000E-05 

23.688700 

0. 27067 300E-04 

‘ 

23,078200 

0.25862200E-04 

22 . 503200 

0.25311400E-04 

40 

21.977400 

0.2412240OE-04 

n 

21.257700 

21.255900 

G.24105100E-04 

0.23480700E-04 

20.763700 

0.22769600E-04 

20.418000 

0.21503500E-04 

19.954900 

0.19712600E-04 

m 

19.455400 

0.18627100E-04 

19.004100 

0.17825600E“04 

18.708400 

0.16949600E-04 

18.616500 

0.15456700E-04 

18.188900 

0.14153100E-04 

,*  *  ‘  i 

17.590300 

0.15377500E-04 

'  ‘ 

17.743100 

0.13075600E-04 

V 

17.661200 

0.12305100E-04 

w... 

17.430100 

0.13105500E-04 

17.173600 

0.13359400E-04 

16.975200 

0.12725900E-04 

",  .\"i 

16.431200 

0.12188600E-04 

16.579500 

0.11153500E-04 

li.rfi* 

o 

16.302900 

0.10068900E-04 

16.168200 

0.93669000E-05 

-V'V. ; 

15.782500 

0.92354000E-05 

15.879800 

0.92480000E-05 

15.754300 

0.87858000E-05 

w 

15.542300 

0.89842000E-05 

15.396700 

0 . 80098000E-05 

15.211100 

0.66224000E-05 

15.101000 

0.57909000E-05 

’  y 

15.153800 

0.49209000E-05 

14.946400 

0.49933000E-05 

15.001900 

0.53894000E-05 

"'-1 

14.728000 

0.45468000E-05 

0  1 

14.842600 

0.8297600QE-03 

15.145600 

0.10003500E-04 

14.879900 

0 . 88311000E-05 

"  /  .  ’■] 

13.872800 

0.78366000E-05 

- 

13.823700 

0.67594000E-05 

; 

13.657200 

0.54952000E-05 

i  ... 

13.174000 

0.49055000E-05 

'  .  '  ,  1 

13.073100 

0.4714200CE-05 

12.982100 

0.45890000E-05 

12.846500 

0 . 49453000E-05 

12.554500 

0.50342000E-05 

v 

t 

12.684600 

0.48464000E-05 

i 

12.543600 

0 . 48563000E-05 
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Specimen  83270  (R“.l, 

4K  de/dK 

14.753500  0.31376710E-03 

14.790800  0.25961760E-03 

14.809900  0.22712950E-03 

14.835400  0.17298000E-03 

14.843600  0.29995220E-03 

14.882700  0.30562930E-03 

14.917200  0.21157040E-03 

14.959100  0.25876720E-03 

14.991800  0.31683790E-03 

15.037300  0.27470810E-03 

15.052800  0.17003120E-03 

15.090100  0.17757440E-03 

15.109200  0.25813200E-03 

15.152900  0.27185380E-03 

15.216600  0.31475920E-03 

15.263900  0.34575130E-03 

15.315700  0. 29609390 E-03 

15.339400  0.17795240B-03 

15.358500  0.1989681OE-O3 

15.382200  0.12462180E-03 

15.390300  0.21940900E-03 

15.407600  0.25294040E-03 

15.457700  0.35510560E-03 

15.517700  0.40681020E-03 

15.581400  0.36412920E-03 

15.643300  0.37520000E-03 

15.666900  0.28497970E-03 

15.704200  0.28662150E-03 

15.734300  0.31705450E-03 

15.786100  0. 307243 50E-03 

15.834300  0.34959380E-03 

15.883500  0. 281822 30E-03 

15.912600  0.30527100E-03 

15.954400  0.27892860E-03 

15.993500  0.19253110K-03 

16.024500  0.18644060E-03 

16.029900  0.18772000E-03 

16.067200  0.29848760E-03 

16.131800  0.40602280E-03 

16.215500  0. 440471 60E-03 

16.279200  0. 38465670E-03 

16.315600  0.38350320E-03 

16.370200  0.34325130E-03 

16.398400  0.33034190E-03 

16.457600  0.34993240E-03 

16.516700  0.33201510E-03 

16.558500  0.30278290E-03 

16.579500  0.29382620E-03 

16.626800  0.36095600E-03 

16.668600  0.37718030E-03 


y-1  Hz ,  t jj-50  tec) 

AK  d»/dN 

16.811500  0.41039290E-03 

16.866100  0.35951500E-03 

16.907900  0.28622780E-03 

16.931600  0.32201120E-03 

16.973400  0.34307810E-03 

17.027100  0.41322750E-03 

17.099000  0.42204640E-03 

17.168100  0.42090470E-03 

17.228200  0.42271570E-03 

17.321900  0.43492040E-03 

17.392800  0.38186540E-03 

17.456500  0.40452680E-03 

17.506600  0.42354250E-03 

17.566600  0.38299530E-03 

17.622100  0.36844810E-03 

17.672100  0.32557810E-03 

17.729500  0.38605040E-03 

17.768600  0.35463710E-03 

17.890500  0.51704620E-03 

17.984200  0.47909350E-03 

18.071600  0.45873920E-03 

18.131600  0.38506620E-03 

18.173500  0.37372370E-03 

18.215300  0.32190490E-03 

18.247100  0.40873930E-03 

18.334500  0.50743990E-03 

18.420900  0. 51023520E-03 

18.504600  0.51822730E-03 

18.583800  0.4IS578650E-03 

18.649300  0.50759740E-03 

18.792100  0.48326680E-03 

18.870400  0.49897540E-03 

18.938600  0.55271540E-03 

19.030500  0.59952640E-03 

19.119700  0.56747920E-03 

19.227000  0.58948700E-03 

19.327100  0.57779410E-03 

19.410800  0.51618010E-03 

19.470800  0.48027460E-03 

19.547300  0.50232180E-03 

19.630100  0.55598310E-03 

19.722900  0.58472320E-03 

19.819300  0.63842390E-03 

19.936700  0. 64791210E-03 

20.041300  0.64499870E-03 

20.151400  0.64145540E-03 

20.251500  0.62602240E-03 

20.353400  0.63366020E-03 

20.453400  0  63783340E-03 

20.569900  0.66936870E-03 


C 


57 


Specimen  83270  (R-.l.V-l  Hz,  tg-50  see) 

At  da/dN 

21,037500  0.72157340E-03 

21.167600  0.72889620E-03 

21.294100  0.75157330E-03 

21.429700  0.76110080E-03 

21.575200  0.78145510E-03 

21.710800  0.75358120E-03 

21.849100  0.73637650E-03 

21.980100  0.76129770E-03 

22.114800  0. 750077 20E-03 

22.254900  0.77858110E-03 

22.398600  0.82381720E-03 

22.563300  0.86696680E-03 

22.738000  0.88318720E-03 

22.905400  0.88822660B-03 

23.090100  0.90176980E-03 

23.264800  0.88283290E-03 

23.428500  0.8P783290E-03 

23.611400  0.92484070E-03 

23.800600  0.97708470E-03 

24.000800  0.10372814E-02 

24.478400  0.11515331E-02 

24.737700  0.11334100E-02 

24.996100  0.11277112E-02 

25.470100  0.11022813E-02 

25.713100  0.11258639E-02 

25.979600  0.11760606E-02 

26.267100  0.12541707B-02 

26.572800  0.13001549E-02 

26.889400  0 .13695642E-02 

27.246100  0 .14414932Er-02 

27.622700  0 .15011781E-02 

28.024900  0.15573197E-02 

28.445200  0.16019653E-02 

28.896500  0.16462172E-02 

29.352300  0.16892486E-02 

29.856300  0.17250753E-02 

30.378600  0.17315713E-02 

30.899000  0.17843271E-02 

31.455800  0.18458625E-02 

32.049000  0.19523583E-02 

32.714000  0.20412558E-02 

33.479200  0.21142871EH)2 

34.277100  0.21683028E-02 
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Specimen  83288 

(First  Proof  Test, 

R-.5,  V-.01 

He,  tg**50  sec) 

• 

4K 

di/dN 

d./dN 

10.0516 

3.88838E-06 

10.9195 

2.27606E-04 

10,0534 

8.14447E-05 

10.9796 

.0005924 

10.0598 

.000159 

11.056 

5.94763E-04 

10.0943 

8.46061E-05 

11.1261 

5.05235E-04 

10.0989 

1.5274E-04 

11.177 

4.52716E-04 

* 

10.0534 

1.63189E-04 

11.1688 

4.14527E-04 

10.0852 

3.U681E-04 

11.2325 

3.71373E-04 

10.1262 

1.86389E-04 

11.2598 

2.13909E-04 

10.1153 

7.42676E-05 

11.2826 

3 .73346E-Q4 

10.1025 

4.34054E-04 

11.3108 

4.80944E-04 

10.1307 

1 .15413E-04 

11.3481 

4.66968E-04 

# 

10.1371 

3.27425E-C4 

11.3826 

5 . 28149E-04 

10.1662 

5.58739E-04 

11.4927 

5.9055E-04 

10.2153 

5.91960E-04 

11.5528 

4.38188E-04 

10.2408 

1.0937E-04 

11.5555 

2.78582E-04 

10.2363 

2.83673E-04 

11.5691 

2.13885E-04 

10.2153 

7.32676E-05 

11.5746 

3.71094E-04 

0 

10.2144 

2.80799E-04 

11.6019 

5.3053E-04 

10.2736 

1.1646E-04 

11.6947 

7.23306E-04 

10.2927 

2.30567E-04 

11.7802 

7.18542E-04 

10.259 

4.02283E-04 

11.8585 

6.11889E-04 

10.2572 

3.53208E-04 

11.9149 

6.15274E-04 

10.2863 

5.93503E-04 

11.9558 

3.73374E-04 

o 

10,3491 

8.74172E-05 

11.9813 

3.12035E-04 

10,3591 

3.64913E-04 

12.0104 

3. 5661 E-04 

10.3727 

3.71133E-04 

12.0386 

4.2992E-04 

10,3318 

1.97468E-04 

12,0795 

5.44605E-04 

10,3682 

3 .42086E-04 

12.1423 

4.74723E-04 

10,4155 

2.48755E-04 

12.2197 

5.27834E-04 

10.4701 

3 .74212E-04 

12.2542 

5.58188E-04 

(D 

10,3764 

2.86594E-04 

12.2897 

6.85471E-04 

10,41 

6.24763E-04 

12.3752 

5.64802E-04 

10.4665 

5.95589E-04 

12.4353 

6.10038E-04 

10.5292 

5.27873E-04 

12,5217 

7.08975E-04 

10.491 

2.01511 E-04 

12.5818 

6.88857E-04 

10.5001 

4.22086E-05 

12.6473 

4.0681E-04 

10.4892 

1.54141E-04 

12.6591 

5.29369E-04 

10.5292 

3 .11598E-04 

12.7464 

7.29684E-04 

' 

10.5793 

4.82007E-04 

12.7992 

7.36888E-04 

10.5975 

3.01027E-04 

12.9011 

6.39841E-04 

10.6639 

4. 2 7401 E-04 

12.9748 

6.60117E-04 

10,6857 

3.39507E-04 

13.0794 

6.09093E-04 

V 

10.7039 

2.74098E-04 

13.0776 

4.23306E-04 

10.7394 

2.14059E-04 

13,1095 

5.8429E-04 

10.7339 

3.12889E-04 

13.1932 

4.66889E-04 

10.7703 

4.77991E-04 

13.265 

6.44959E-04 

10.8158 

5.98424E-04 

13.3151 

6.69684E-04 

10.9004 

4.31023E-04 

13.4106 

7.14408E-04 

1 

10.9578 

3.00059E-04 

13,4725 

5.82164E-04 

10.9796 

1.96527E-04 

13,518 

7.58818E-04 

59 


5. 


Specimen  83288 

(First  Proof  Test, 

K- . 5 , V  -.01 

Hz,  tg-50  sec) 

AK 

ds/dN 

AE 

ds/dN 

13.7181 

6.11062E-04 

25.0107 

2.60377E-03 

13.7854 

5.43188E-04 

25.4492 

.0027085 

13.8546 

5.4681E-04 

25.9123 

2.85177E-03 

13.8928 

6.00904E-04 

26.4582 

3.10P84E-03 

13.9237 

6.6807E-04 

27.1114 

3.2962iE-03 

14.0593 

8.970B5E-04 

27.8338 

3.23137E-03 

14.1612 

9.20471E-04 

28.4989 

3.10948E-03 

14.2913 

1.03264E-03 

29.1712 

3.14496E-03 

14.4278 

9.42833E- 04 

29.8345 

3.12881E-03 

14.5288 

7.67951E-04 

30.5214 

3.19814E-03 

14.6061 

7.06967E-04 

31.3302 

3.44704E-03 

14.688 

7.69526E-04 

32.2737 

3.86003E-03 

14.7626 

7.60983E-04 

33.3518 

4.12716E-03 

14.8599 

7.72636E-04 

34.5646 

4.44645E-03 

14.9837 

8.52518E-04 

36.0212 

4.76298E-03 

15.101 

8.38896E-04 

37.6252 

5.17834E-03 

15.1838 

8.8492E-04 

39.5422 

5.49645E-03 

15.2839 

9.45982E-04 

15.4222 

1.01846E-03 

15.5841 

1.04437E-03 

15.7379 

9.94329E-04 

15.8671 

9.37715E-04 

15.9672 

7.90628E-04 

16.0518 

8.89408E-04 

16.1755 

1.02583E-03 

16.3402 

1.13244».-03 

16.5203 

1.23267E-03 

16.7132 

1.29996E-03 

16.9134 

1.41551E-03 

17.1172 

1.33114E-03 

17.3373 

.0013359 

17.563 

1.30027E-03 

17.7568 

1.24279E-03 

17.9487 

1.15331E-03 

18.0943 

1.23291E-03 

18.329 

1.32106E-03 

18.561 

1.38665E-03 

18.8067 

1.56035E-03 

19.0978 

1.68019E-03 

19.4381 

1.76267E-03 

19.7583 

1.75098E-03 

20.1468 

1.83937E-03 

20.5171 

1.76511E-03 

20.8383 

1.74874E-03 

21.2104 

1 . 805S2E-03 

21.6371 

1.86945E-03 

22.0338 

1.97078E-03 

22.5023 

2.14645E-03 

23.0591 

2.28842E-03 

23.646 

2.34996E-03 

60 


Specimen  83289  (Second  Proof 

AK  da/dN 

9.5666522  0.16375164E-03 

9.6330684  0.12005881E-03 

9.6785589  0.76362050E-04 

9.7040335  0.32670406E-04 

9.6794687  0.28036164E-05 

9.6794687  0.29657420E-05 

9. 67855^9  0 . 38258190E-05 

9.6894766  0.58436889E-05 

9.6894766  0.83944712E-05 

9.69 49354  0.52149501E-05 

9.6931158  0.83649437E-06 

9.6894766  0.16207054E-04 

9.7049433  0.99232082E-05 

9.7122218  0. 149247 73E-04 

9.7140414  0.16815320E-04 

9.7276886  0.19088938E-04 

9.7331475  0.18692482E-04 

9.7386063  0.11940133E-04 

9.7595319  0.15075166E-04 

9.7477044  0.37350712E-05 

9.7504338  0.87334469E-05 

9.7568025  0.95180909E-05 

9.7622614  0. 1993381 8E-04 

9.7777281  0.93503747E-05 

9.7731791  0.17951145E-05 

9.7786380  0.84440774E-05 

9.7868262  0.14115719E-04 

9.7813674  0.19131457E-05 

9.7850066  0 .24303100E-05 

9.7950145  0.12645250E-04 

9.7968342  0 .43988100E-05 

9.7950145  0.61940819E-05 

9.7977440  0 . 60645546E-05 

9.8050224  0 .71054974E-05 

9,8095715  0.52429028E-05 

9.7804576  0.61763654E-05 

9.7813674  0.14226349E-04 

9.7804576  0.U822023E-04 

9.7968342  0.15748787E-04 

9.8050224  0.24609793E-04 

9.8195794  0 . 17135005E-04 

9.8032028  C .77669134E-05 

9.8032028  0. S8507674E-05 

9.8268579  0 .20279486E-04 

9.8496031  0.20324368E-04 

9.8514228  0.50755803E-05 

11.060560  0.21803105E-04 

11.062380  0.80696687E-04 

11.096953  0.889O1395E-04 

11.140624  0.68696711E-04 


Test,  R**.25,V»1  Hz,  tg-10  sec) 

AK  da/dN 

11.550948  0.49716435E-04 

11.574603  0.64177035E-04 

11.610995  0.66799077E-04 

11.662854  0.77932913E-04 

11.716533  0.77503780E-04 

11.753835  0.67680965E-04 

11.785679  0.65306954E-04 

11.814793  0.58877833E-04 

11.849366  0.69960488E-04 

11.893946  0.66799077E-04 

11.940347  0.55724296E-04 

11.963092  0.56110122E-04 

11.993116  0.43354243E-04 

12.009492  0.47873919E-04 

12.028598  0.497 75490E-04 

12.067720  0.603 89641E-04 

12.103203  0.80110074E-04 

12.146873  0.77850236E-04 

12.209650  0.83751799E-04 

12.254231  0.70212456E-04 

12.296082  0.64133728E-04 

12.325196  0.6141 7198E-04 

12.342483  0.57574686E-04 

9.5757503  0.77125828E-04 

12.704587  0.84940773E-04 

12.752807  0.77669134E-04 

12.797388  0.85090379E-04 

12.842878  0.91799027E-04 

12.909294  0.91113989E-04 

12.958424  0.87724232E-04 

13.015742  0.91657295E-04 

13.075789  0.86554943E-04 

13.114911  0 . 79763618E-04 

13.156763  0.86110062E-04 

13.217720  0. 97936809E-04 

13.276857  0.96413190E-04 

13.342364  0.95125791E-04 

13.406050  0.97696652E-04 

13.459729  0.91995877E-04 

13.513408  0 . 82110070E-04 

13.560718  0.83984082E-04 

13.613487  0.87476201E-04 

13.667166  0 . 94062S01E-04 

13.722664  0.98456493E-04 

14.012894  0.81507709E-04 

14.069302  0.39 936927E-04 

13.725394  0.55417210E-03 

13.774524  0.28384588E-03 

13.931011  0.25520814E-03 

14.081129  0.22632631E-03 
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Spocimen  83289 
AX 

14.359531 

14.449503 

14.510560 

14.570607 

14.647031 

14.742561 

14.842640 

14.925433 

15.012775 

15.092838 

15.184729 

15.286628 

15.382158 

15.484966 

14.448693 

14.414120 

14.463250 

14.556960 

14.637023 

14.723455 

14.821715 

14.901778 

14.978202 

15.071913 

15.143337 

15.240227 

15.344856 

15.448574 

15.556841 

15.666019 

15.764278 

15.867087 

15.973534 

16.088170 

16.200077 

16.322901 

16.452094 

16.534887 

16.624048 

16.729586 

16.796003 

16.858779 

16.919737 

16.978874 

17.034373 

17.111707 

17.182672 

17.253637 

17.321873 

17.381010 


(Seoond  Proof 
da/dN 

0.16236975E-03 

0.13773988E-03 

0.11799976E-03 

0.11170450E-03 

0.12144464E-03 

0.13218871E-03 

0.13134225E-03 

0.12788950E-03 

0.12346432E-03 

0.125291082-03 

0.12902336E-03 

0.13247611E-03 

0.13970050E-03 

0.14270444E  -03 

0.22709403E-04 

0.26082231E-04 

0.86078566E-04 

0.11127930E-03 

0 .12012574E-03 

0.12584620E-03 

0.12479502E-03 

0.11966511E-03 

0 . 11915330E-03 

0.12270841E-03 

0.12481864E-03 

0.13082257E-03 

0.13948791E-03 

0.13904303E-03 

0.13973200E-C3 

0.14017294E-03 

0.13651941E-03 

0.13977924E-03 

0.14208239E-03 

0.14489341E-03 

0.15164930E-03 

0.16039338E-05 

0.16968863E-03 

0.17641697E-03 

0.19575945E-03 

0 .18638151E-03 

0.17344453E-03 

0.17692090E-03 

0 . 17616893E-03 

0.17162170E-03 

0.17344059E-03 

0.19096418E-03 

0.19263740E-03 

0.17751145E-03 

0.17065320E-03 

0.17100359E-03 


Test,  E-.25.V-1 
AK 

17.503835 

17.564792 

17.627569 

17.703083 

17.783146 

17.877767 

17.962379 

18.033.744 

18.105219 

18.186192 

18.278993 

18.371793 

18.460045 

18.549206 

18.646556 

18.731168 

18.825789 

18.914040 

19.022308 

19.137853 

19,237933 

19.338922 

19.439001 

19.536350 

19.653716 

<9.771081 

19.872070 

19.98670 6 

20.096793 

20.216888 

20.345171 

20.468906 

20.587181 

20.717284 

20.841928 

20.990227 

21.150353 

21.317758 

21.477885 

21.643470 

21.800867 

21.961904 

22.141136 

22.325828 

22.511429 

22.710677 

22.917204 

23.112813 

23.302964 

23.481286 


He,  tg-10  seo) 

<W  dN 

0.16316896E-03 

0.16439730E-03 

0 .192247 64E-03 

0.20194447 £-03 

0.20551927E-03 

Q.20424368E-03 

0.20477518E-03 

0.20538147E-03 

0.19379095E-03 

0.20825548E-03 

0.21360587E-03 

0.22461372E-03 

0. 2200625 5E-0 3 

0.21486571E-03 

0.21633814E-03 

0.21813342E-03 

0.22385781E-03 

0.23236173E-03 

0.23813731E-03 

0.23318850E-03 

0. 23875 542E-03 

0.22896410E-03 

0.23054284E-03 

0.23260189E-03 

0.24056250E-03 

0.24128298E-03 

0.23469244E-03 

0.24292864E-03 

0.24705462E-03 

0.24826721E-03 

0.25256642E-03 

0.25299555E-03 

0.25893648E-03 

0. 2601569 5E-03 

0.27507425E-03 

0.29486948E-03 

0.30531040E-03 

0.30379860E-03 

0.30285765E-03 

0.30936551E-03 

0.29999152E-03 

0.31017260E-03 

0.31699148E-03 

0.32900327E-03 

0. 33371586E-03 

0.33447176E-03 

0.33157019E-03 

0.31678676E-03 

0 . 31501511E-03 

0.31123559E-03 
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Specimen  83289  (Second  Proof  Test,  R-.25,  V»1  Hz,  tg-10  sec) 


At 

23.87(144 

24.095408 

24.334688 

24.585796 

24.862378 

25.152607 

25.457394 

25.743074 

26.048770 

26.370843 

26.702924 

27.052291 

27.405297 

27.748295 

28,108580 

28,501618 

28.916491 

29.377765 

29.881800 

30.391293 

30.908975 

31.433936 

32.015304 

32.670368 

33.386388 

34.120605 

34.822978 

35.481680 

36.189513 

37.021079 

38.016411 

39.176419 

40,430137 

41.779385 

43.244179 

44.847265 

46.312454 


da/dN 

0.33266468E-03 
0.3542 7881E-03 
0.37S65672E-03 
0.39255432E-03 
C.40555036E-03 
0.41334562E-03 
0.416770S1E-03 
0.42Q94408E-03 
0.42637709E-03 
0.4271 644  9E--03 
0.43712510E-03 
0,4392117 IE-03 
0.43925108E-03 
0.43736132E-03 
0.45495971E-03 
0.47499904E-03 
0.49358168E-03 
0.51964462E-03 
0.52881783E-03 
0.53618002E-03 
0.53964458E-03 
0.56436893E-03 
0.58826652E-03 
0.60369956E-03 
0.61684914E-03 
0.60216413E-03 
0.59188856E-03 
0.58995943E-03 
0.62716408E-03 
0.68621908E-03 
0.74106149E-03 
0.79708500E-03 
0.83464398E-03 
0 . 86613998E-03 
0.88704545E-03 
0.91527374E-03 
0.84704553E-03 


NOTE:  All  of  above  date  (83264P,  83266G,  83270B,  _83270,  83288, 

83289)  is  in  English  Units.  *  Ksi  '/inches, 
da/dN  *  inches  /cycle. 


and 


Constant  Ai l 

IN  718  da/dN  Data,  Variable  Bold  Time 


(T-1200F,  B-0.1,  V-lHz) 


AK-22.745  Kai  VIT 


Km  -25.3  Kai  Via" 

Dll 


Hold  Time  (») 


0.0 

10.0 

25.0 

50.0 

200.0 

500.0 


1.55  E-5 
7.48  E-5 
1.61  E-4 
2.80  E-4 
1.02  E-3 
2.80  E-3 


4K-32.75  Kai  Via 


-36.4  Kai  Via" 


0.0 

1.0 

4.0 

10.0 

30.0 

50.0 

100.0 

200.0 

300.0 

500.0 


3.78  E-5 
5.08  E-5 
8.86  E-5 

1.93  E-4 
4.76  E-4 
7.18  E-4 
1.10  E-3 
2.14  E-3 

3.94  E-3 
6.65  E-3 


t?*-  -"n  «Wv  *,"■ 


n  n «  non  nno  <100  000  rjooonoooo 


Appendix  B:  FORTRAN  Code  For  SINH  Interpoletive  Model 


PROGRAM  SINB718 

THIS  PROGRAM  GENERATES  THE  CONSTANTS  FOR  THE  HYPERBOLIC 
SINE  EQUATION  AND  COMPUTES  CRACK  GROWTH  RATES  FOR 
DESIRED  STRESS  INTENSITY  RANGES.  THE  MODEL  IS  ISOTHERMAL 
(1200F) ,  AND  IS  SPECIFIC  TO  INCONEL  718.  REQUIRED  INPUTS 
ARE  STRESS  RATIO.  FREQUENCY  (HZ).  AND  HOLD  TIME  (SEC). 

READ  INPUT 

PRINT  * . ' INPUT  STRESS  RATIO;' 

READ  *,R 

PRINT  *. 'INPUT  FREQUENCY  (HZ);' 

READ  +.FBEQ 

PRINT  *» 'INPUT  HOLD  TIME  (SEC);' 

READ  *, THOLD 

SINH  CONSTANT  Cl 

C1-.5 

SINH  CONSTANT  C2 

C2-3 . 1 97-1 . 1675*AL0G10 ( ( 1 . -R> / . 9 ) - . 1505*ALOG10 ( FREQ) 

•  +.2O53*ALOG10(THOLD+l.) 

SINH  CONSTANT  C3 

C3— 1 , 721- . 6243*ALOG10 ( ( 1 .-R) / . 9 ) +  .00717  *ALOG10 ( FREQ) 

•  + ,0677*ALOG10 (THOLD+1 . ) ) 

SINH  CONSTANT  C4 

C4— 3 .935+1 .0131*ALOG10(  (1  .-R)  / ,  9)  -  ,523*ALOG10  (FREQ) 

•  + . 8127  *ALOG10 ( THOLD+1 . ) ) 

OUTPUT 

PRINT  20.C1.C2.C3.C4 
20  F0RMAT(//2X,'C1«',F10.4,/2X, 'C2-' .F10.4./2X. 'C3-',F10.4, 

•  /2X, *C4“' .Flo. 4) 

PRINT  R  FREQ  THOLD' 

PRINT  30. R, FREQ, THOLD 
30  F0RMAT(2X.3F10,4, //) 

PRINT  •. 'DESIRE  DA/DN  VALUES 7  1  FOR  YES,  0  FOR  NO' 

READ  + , IFLAG 
IF(IFLAG.EQ.O)  GOTO  50 

40  PRINT  *, 'INPUT  DESIRED  DELTA  K  (KSI  ROOT  IN)’ 

READ  +.DK 
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DADN-10**<C1*SINH(C2*(AL0G10(DK)+C3))+C4) 

PRINT  42.DADN.dk 

FORMAT(2X, 'CRACK  GROWTH  RATE  IN  IN/ CYCLE  IS: ' , 
£12.5,31, 'FOR  DI«',F6.2,/) 

PRINT  •.'ANOTHER  DA/DN  VALUE?  1  FOR  YES,  0  FOR  NO' 

READ  • , JFLAG 

IF( JFLAG.EQ.l)  GOTO  40 

PRINT  •, 'CHANGE  TEST  PARAMETERS?  1  FOR  YES,  0  FOR  NO 
READ  * ,  KFLAG 
IF (KFLAG. EO.l)  GOTO  10 


ouooyuuuu  you  u  u  u  u  u  u  u  u  u  u  u  u  u  u  u  wow 


Appendix  C:  FORTRAN  Code  For  USE  Interpolative  Model 


PROGRAM  MSE718 

THIS  PROGRAM  GENERATES  THE  CONSTANTS  FOR  THE  MODIFIED 
SIGMOIDAL  EQUATION,  AND  COMPUTES  CRACK  GROWTH  RATES 
FOR  DESIRED  STRESS  INTENSITY  RANGES.  THE  MODEL  IS  ISOTHERMAL 
(1200F),  AND  IS  SPECIFIC  TO  INCONEL  718.  REQUIRED  INPUTS  ARE 
STRESS  RATIO,  FREQUENCY  (HZ),  AND  HOLD  TIME  (SEC). 

READ  INPUT 

PRINT  *, 'INPUT  STRESS  RATIO:1 
READ  *  & 

PRINT  *, 'INPUT  FREQUENCY  (HZ):' 

READ  * , FREQ 

PRINT  *, 'INPUT  HOLD  TIME  (SEC):' 

READ  * .THOLD 

UPPER  ASYMTOTE  (ASSUMED  CONSTANT) 

DKC-109.9 
LOWER  ASYMTOTE 

DKSTAR-10** (1 . 07918- .048646*ALOG10 (FREQ) 

*  +.05675*ALOG10(THOLD+1.) ) 

VERTICAL  LOCATION  OF  INFLECTION  POINT 

DADNI-10 . 0** (-3 . 933+1 . 0131*AL0610 ( (1 .-R) / . 9 ) - . 523*AL0G10 (FREQ) 

*  +. 8127*ALOG10(THOLD+1 . ) ) 

HORIZONTAL  LOCATION  OF  INFLECTION  POINT 

DKI»10.0'**(1 ,721+.6243*ALOG10( (1 .-R) / .9)-. 00717 *ALOG10( FREQ) 

*  -.0677*AL0G10(THOLD+l . ) ) 

SLOPE  OF  INFLECTION, POINT 

DADNIP-1 . 7-1 . 8451*ALOtlO ( ( 1 . -R) / . 9 ) - . 2 18 5*ALOG10 ( FREQ) 

*  +.1956*ALOG10(THOLD+1 , ) 

LOWER  SHAPING  COEFFICIENT  (ASSUMED  CONSTANT) 

Q-1.0 

UPPER  SHAPING  COEFFICIENT 

(.  D**-(SQRT(Q)  *ALOG(DKI/DKC)  /ALOG(DKI/DKSTAR) )  **2 

C 
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OTHER  COEFFICIENTS  8'  AMD  P 


CONST-ALOG ( DEI/DKSTAR) 

CONST1-AIQG ( DEC/DEI) 

BPRIME-ALOG (DADNI ) -Q*ALOG( CONST) -DIALOG ( CONST1 ) 
P“DADNIP~Q/ CONST+D/ CONST! 

OUTPUT 

PRINT  20 , BPRINE, P . 0. D , DESTAR , DEC , DADNI , DADNIP 
FORMAT(//2X, 'B' .F10.4./2X, 'P-* ,F10.4»/2X, 'Q-' ,F10.4, 
/2X, 'D"' ,F10.4,/2X, 'DELTA  K*-' ,F10,4,/2X, ’DELTA  EC-*, 
F10.4, /2X, 'DA/DNI-' ,E15.5,/2X, 'DA/DNI* *«» ,F10.4,//) 
PRINT*.*  R  FREQ  THOLD* 

PRINT  30, R, FREQ, THOLD 
FOPMAT( 2X, 3F10 ,  4 ,  /  / ) 

PRINT  *, 'DESIRE  DA/DN  VALUES/  1  FOR  YES,  0  FOR  NO' 

READ  *,IFLAG 
IF(IFLAG.EQ.O)  GOTO  30 

PRINT  * , *  INPUT  DESIRED  DELTA  E  (ESI  ROOT  IN)' 

READ  *,DE 

DADN-EXP(BPRIME) * ( (DE/DEI) **P) * ( ALOG ( DE/ DESTAR ) **Q) 

* ( ALOG ( DEC/DK) **D ) 

PRINT  42.DADN.DE 

FORUAT(2X, 'CRACK  GROWTH  HATE  IN  IN/ CYCLE  IS: » , 

E12.5.3X, 'FOR  DE-',F«.2,/) 

PRINT  *, 'ANOTHER  DA/DN  VALUE?  1  FOR  YES,  0  FOR  NO* 

READ  *,JFLAG 
IF(JFLAG.EQ.l)  GOTO  40 

PRINT  * , *  CHANGE  TEST  PARAMETERS?  1  FOR  YES,  0  FOR  NO' 
READ  *,EFLAG 
IFdFLAG.EQ.l)  GOTO  10 
STOP 
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Thi»,,^wrk-p  re  tents  in  examination  of  interpolates  schemes  for 
fatigue  crack  growth  ut  elevated  temperature  (1200 *F).  An  interpolative 
scheme  involving  the  linear  superposition  of  effects  due  to  stress 
ratio,  loading  frequency,  and  peak  load  dwell  ox  hold  time,  is  applied 
to  two  state-of-the-art  crack  growth  rate  prcdlotion  models^  These 
models  arc  the  SINE  model,  developed  by  Pratt  an^iThltnieyA^roraft,  and 
the  MSE  model,  developed  by  General  Eloo  tr  io.1*  SINt^fsbisHed  on  the 
hyperbolic  sine  function,  and  MSfylsoased  on  a  sigmoidal  equation. 

The  results  of  an  experimental  program  are  presented.  Fatigue 
crack  growth  rate  teats  are  performed  on  compact  tension  specimens 
according  to  A3TM  standard!.  Two  interpolative  models  are  developed 
from  the  resulting  data  base.  Additional  teats  are  performed  and  the 
model  predictions  are  compared  to  these  additional  tests. 

It  is  found  that  linear  (on  log-log  aoale)  functional  forms  can  be 
assumed  to  relate  model  constants  to  the  teat  parameters  (load  ratio, 
frequency,  and  hold  time).  Also,  it  is  found  that  these  functional 
forms  can  be  applied  to  both  mode  la.  v  More  work  is  needed  to  determine 
whether  the  effects  of  variations  in  eaoh  uf  the  test  parameters  on 
crack  growth  rate  are  independent  of  eaoh  other  and  can  be  linearly 
superimposed. 

When  the  same  functional  forma  are  applied  to  tho  SINH  and  the  MSE 
models,  they  behave  much  alike.  '  C**f. 
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